UNIVERSITY OF ILLINOIS 
DIGITAL COMPUTER 



ILLIHOIS CODE M 10 - 153 



TITLE 

TYPE 

METHOD OF USE 



RmCHUfG OF TAJPES 



Estimation of Communalities by the Method of 
Maximum Likelihood. 
Entire program. 

The program is read into the memory in the usual way. 
After a short time, Illiac will stop. Bead in para- 
meters and the elements of the correlation matrix. 
•Hie computer will stop and is now prepared to read in 
the estimated elements of the matrix 0* Beplace the 
original tape in the reader and restart the computer. 
The "black switch Bust be ia the ignore position* If 
convergence has not occurred in the allotted time* 
raise the black switch to the obey position. After a 
short time the computer will come to a stop. At that 
ti»e raise the white switch up and down. This will 
read in the print routine and the results will be 
printed out shortly* 

At the beginning of each data tape four parameters are 
necessary, ^hey are as follows j 

(1) Let s = size of sample. Then place sS on the tape- 

(2 ) Let f = the number of estimated factors . Place 
fF on the tape • ' 

(3) Let n = the number of variables. Place nft on the 
'■tape. 

(k) Let j = the number of places to be printed. Place 

jj on the tape. 
Since the chi-square value varies with the number of 
factors t it is necessary to provide several values. 
These can be found in Table I for the 5 ^ level. ' 



-2- 

Each value (which must be sealed "by 10 ) is punched 
as a sign followed by up to 12 decimal digits* The 

final value in the sequence is terminated by an N. 

-2 
The values of the correlation matrix scaled "by 10 

are next punched on the tape. The lower off -diagonal 

elements and diagonal elements are punched, row "by row, 

as a sign followed by up to 12 decimal digits • The 

final fraction in the sequence is terminated by an N. 

It is necessary to make some estimation of the matrix G« 

The values scaled by 10"" 1 * should be punehed as a sign 

followed by up to 12 decimal digits* The last value 

is followed by an N, 

MATHEMATICAL The essential steps in the iterative process of estimation 

METHOD and tests of significance may be summarised as follows :. 

(1) Start with a trial dlagoml mtrtx G and tibtala 

the solution to the determin&ntal equation |<2RGr - XI | - 0. 

*0iis glv$s the first approximation to the diagonal elements 

of G. 

g ± */l +^6 U 2 (\/\ - 1 ) + -4[2 2 CV X e ~ X>) + "* + 

where £„,«**, j%L are the f (number of factors assumed) latent 
normalized vectors corresponding to the f largest roots 



• 9 \„ and 



fe f+l + Vi-2 + --«+\)/a-f 



(2) This process is repeated till g. converges to about 
four decimal places. When g converges the sum 
2 



^TIM&TIGN OF G 



TEE ERIHT-OUT 



(3) The estimated factor loading's at any stage are 



^2 

The estimated commonalities are, at any stage, 
2 



h. 
1 



(Sj -i)/gi 



(^) The test for a, specified number k of factors is f 
J\^ = - (S-l - (2h + 5)/6 - (2f)/3 )/^n\f +1 + ... + 

•£, n?^ - (n-f ) ^c n X e Z 

which is a chi-sqtuared test on [(n-f )-n»f ]/2 degrees of 
freedom* This test is made whenTV coaverges to about five 
places. If /} is greater than the chi -square value that 
the number of factors in increased "by one. 
Because the convergence of the iterative process is some- 
what slow, it is desirable to choose values of G which will 
minimise the number of iterations* It has been shown* that 
the squared multiple correlations are the minimum estimates 
of sommunallties , Consequently, the st^ma^e root of diagonal 
elements of the inverse of the correlation matrix foam a 
minimum estimates for the matrix G. 

After the first iteration, the latent roots are printed out. 
*» A t -lue, are print* out after each iteration, ^is 
is followed by the elements of the matrix G s the carnmunalities, 
the largest latent roots 3 and the factor loadings* If too 
many factors are being extracted or if there are too few values 
of cM*s|^6 ; then ten F's will "be puncked out and the 
computer will stop. 



*L. Guttman, Multiple Rectilinear Prediction and the Resolution into Components, 
Psychometrika, Vol. 5 (19^0) pp, 75~90 
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!EABLE I 


























decked op Tsurnxm 












f /n 


5 


6 


7 


8 


9 


10 


11 


12 




14 


15 


16 


17 


1 


5 


9 


14 


20 


27 


35 


44 


54 


65 


77 


90 


io4 


119 


2 


1 


4 


8 


13 


19 


26 


34 


43 


53 


64 


76 


89 


103 


3 






.3 


7 


12 


l8 


25 


33 


42 


52 


63 


75 


88 


4 








2 


6 


11 


17 


24 


32- 


4i 


51 


62 


74 


5 










1 


5 


10 


16 


23 


31 


40 


50 


61 


6 














4 


9 


15 


22 


30 


39 


49 


T 
















3 . 


8 


14 


21 


29 


38 


8 


















2 


7 


13 


20 


28 


9 




















1 


6 


12 


19 


10 
























5 


11 


11 


























4 



CHI-SQUARE* A3? JJgftKVEL 

1 11.1 • 16,9 23.7 31*4 40,1 49*8 60.2 71.9 84.5 98*2 112,9 128.5 145.2 

2 3.84 9 ? 49 15.5 22.4 30 .1 38.9 48,6 59.0 70.7 83»4 97»l HI .7 127-4 

3 7.81"' 14.1 21.0 28*9 37*7 47*4 57*8 69.5 82.2 95«9 U-0„6 

4 5.«99 12*6 19.7 27.6 36,4 46,2 56.6 68.4 81.1 94.8 

5 3*84: u.l 18,3 26.3 35*2 45*0 55*5 67.2 79.9 

6 9.49 16*9 25.O 33.9 43.8 54.3 66.1 

7 7»8i' 15.5 23.7 32.7 42.6 53.I 

8 5;. 99 14.1 22.4 31.4 41.3 

9 3.84 12.6 21*0 30.1 

10 u.i 19.7 

11 . 9.49 
*When degrees of freeckw axe greater than 30 an approximation of ehi-sqiiare is made. 



LOCATION 


ORDER 




* icqsES page i 




CODE XI 




Decimal Order Input 




00 80K 








CQJUct Jfc_ 




Input a Sequence of Beeiaal Fractions 




00 20K 









41 5F 








kl F 


From 8 




1 


81 4f 

LO 3TL 


Fro® 5 




2 


32 5L 








L4 37L 




lead in parameters and convert 


3 


50 F 
' T^ 37L 






k 


S5 F 
kO F 






5 


26 1L 








k-2 6L 


From 2 




6 


L5 F 








40 ( )F 


By 5 




T 


F5 5F 
ho 5F 






8 


LO 3% 
32 L 






9 


L5 2F 
00 IF 






10 


kO F 
L5 4F 






11 


00 2F 








L4 f 






12 


F4 3TL 
10 39? 






13 


66 3&L 
SI F 




i 


lif 
1 1*5 


IA IF 
40 9F 
L*5 2F 




r 2n+5 2f "? 
Compute f S - 1 - '6 - 3 \ 
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LOCATION 


ORDER 


NOTES PAGE 2 




50 2F 






l6 


7^ 2F 
S5 F 






17 


10 IF 




n(n+l) 




ho 11F 




Compute 2 


18 


IA- 36L 
40 7F 




n(n+l ) 

331 + 2' 


19 


50 2F 
fk 2F 






20 


S5 F 








^0 6f 




2 
el 4- n 


21 


00 20F 
k6 30L 






22 


L5 6F 
lA- 11F 






23 


^0 8F 

00 20F 




3/2 (n 2 + a) 


2k . 


k6 32L 
L5 2F 






25 


LO ij-F 








%0 10F 




n-f 


26 


L5 2F 

46 4f 






2? 


Ik 2F 








40 5F 




2n 


28 


50 988F 








50 28L 




Read In chi~sc[uare values 


29 


26 80F 
22 30L 






30 


50 ( )F 


By 21 






50 30L 




Read in correlation matrix R 


31 


.26 80F 








2^ 32L 




Stop, prepare to read in Matrix G 


32 


50 ( )F 


By 2^ 






50 32L 
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LOCATION 



33 
34 

35 
36 

37 
38 

39 
40 



ORDER 



26 80F 

50 40L 

75 10F 

S5 F 

40 12F 

24 999F 

00 F 

00 331F 

00 F 

00 10F 

00 F 

00 6F 

80 F 

00 4F 

00 F 00 0100 

0000 0000 J 

24 20N 

00 959K 

Code PI 



NOTES 



00 


20K 


Code MO 


00 


152K 


Code Rl 


00 


286K 


Code SI 


00 


9^5K 


92 


131F 


L5 I3F 


52 


74f 


50 


1L 


26 


959F 


L5 


191F 


22 


217F 



.PAfrFi 3 



M 10 



n-f/100 



Print one Nuriber Fractional or 
Integer in a Manner determined by a 
Program parameter 

Eigenvalue - Eigenvector Program 

Square Root Routine 

Natural Logarithm 

Print out chi-square criterion 



LOCATION 



ORDER 
00 F 
00 950K 
kl k$ 
L3 127F 
36 170F 
kl 127F 
kl 5F 
92 lV/F 
92 131F 
L5 (331 )F 
52 72F 
50 kL 
26 9L 
F5 3L 
k2 3L 
F5 5F 
il-0 5F 
LO 278F 
36 3L 
26 170F 

00 161K 
L5 1^7F 
k2 3L 
1^ 3L 
k2 10JL 
kl F 
50 F 
L5 ( )F 
lj-0 ( )F 
F5 F 
40 F 
Kk F 
00 2QF 
Yh 3L 



Fro® 169F 



BOTES 



PAGE k 



M 10 



From 8 
By 6 



Is this first iteration? 

If so, print out all latent roots . 



From 115 



By 1,6 
From 8 
By 0, 6 



Store latent roots in consecutive 
locations • 



LOCATION 


ORDER 
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40 3L 








7 


L5 F 
LO 117L 








3 


36 3L 
26 950F 








9 


4l 5F 
4l 6F 


From 958F 






10 


^9 8F 
4i ?F 


From 30 






XI 


L5 147F 
4a 12L 








12 


46 13L 










LT ( )F 


By 11, 16 


From 18 




13 


L2 ( )F 


By 12, 15 


Find smallest latent roots * 






32 15L 








14 


L5 12L 
00 20F 








15 


46 13L 
F5 12L 


From 13 






16 


42 12L 
F5 7F 








17 


40 7F 
LO 11!?L 








18 


32 12L 
L5 13L 








19 


46 21L 










46 25L 








20 


46 I5F 
41 9F 








21 


L5 ( )F 
40 2F 


From 19 






22 


40 16F 
L4 5F 


> 


- ....... ' .... . ,. . . > 




23 


40 5F 
50 23L 




Sum smallest latent roots 





M 10 



I_ j 






LOCATION 


ORDER 


NOTES PAGE 6 


2k 


22 286F 




Find logarithm of smallest latent 


. 


L5 1^2F 




roots. 


25 


1*0 { )T 


From 19 


-59 
Store 1 - 2 ^ in locations of 




kl 10F 




smallest latent roots. 


26 


50 F 
7J ll8L 


■ 




27 


IA- 120L 

iA 6f 






28 


kO 6F 




Sim logarithms of smallest latent 




F5 ^F 




roots. 


29 


Ao kf 








LO 119L 






30 


32 10L 




■ 




50 S8 




Constant 


31 


L5 5? 
50 31L 






32 


26 286F 








40 7F 






33 


L5 ll6L 








50 33L; 






3^ 


26 286F 
LO 7F 






35 


kO F 

50 F 


; 




36 


7J 118L 








lj-0 F 




• Compute ehi-sqmre criterion 


37 


50 119L 
JO lk2T 






38 


75 P 
00 39F 






39 


Ik 6J 






. 


kO F 






1*0 


50 F 
71 123L 




■ 


in 


LI 13F 
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LOCATION 



WOIES 



PAGE 7 



M 10 





Sh- F 


k-2 


ho F 




S5 F 


h3 


^-0 13F 




19 15F 


kk 


L2 F 




32 if5L 


h5 


26 56L 




L3 (988 )F 


k6 


36 121& 




40 F 


hi 


L7 F 




LO 15F 


k& 


36 56L 




L5 11% 


^9 


LO 145F 




k2 119L 


50 


L5 116L 




LO 121L 


51 


ko 116L 




F5 ^5L 


52 


42 ^5L 




L5 15F 


53 


k6 5^L 




L5 16F 


5^ 


kO ( )F 




L5 5F 


55 


LO 16F 




ho 5F 


% 


26 9A5F 




te 77L 


57 


k-2 8lL 




L5 68L 


58 


22 59L 




L5 3F 



From kk 



By 53 



Has ehi-Square "value converged? 

If not , transfer control to 56 

By 52 

Are too many factors feeing extracted 

If so, transfer to 12^ « 



Does criterion exceed cki -square 
•value? If so, increase number of 
factors "by one, Otherwise, transfer 
control to 56 « 



From ^5, 
From 9^8F 



I4-8 



From 85 



59 


iA 117L 




k2 W 


6o 


kl kF 




L5 1A7F 


61 


k2 62L 




L5 121L 


62 


kO IF 




F7 ( )F 


63 . 


36 64l 




22 73L 


6k 


L5 62L 




te 66l 


65 


iA 3F 




k2 69L 


66 


^2 71L 




50 F 


67 


75 116L 




66 5F 


68 


LI 121L 




Sk SS 


69 


ll-O 2F 




50 ( )F 


TO 


7 J 2F 




AO 2F 


71 


Al F 




50 ( )F 


72 


7J 2F 




iA IF 


73 


AO IF 




F5 62L 


74 


k-2 62L 




F5 ^F 


75 


kO kF 




LO 117L 


76 


32 62L 



HOEES 



PAGE 8 



1 ^ 



By 61, 



From 63 



From 76 



By 65 



By 66 



•Estimate g/ 



From 63 



uoQAnm 



owm 



PAGE 



M 10 





50 76L 


77 


22 152F 




LO ( )F 


78 


40 IF 




19 13F 


19 


L2 IF 




36 8lL 


80 


4l 8f 




36 8lL 


8l 


L5 2F 




40 ( )F 


82 


F5 77L 




42 77L 


83 


42 8lL 




F5 9F 


84 


40 9F 




LO 117L 


85 


32 58L 




L3 &* 


86 


34 87L 




26 999F 


87 


L5 151F 




42 89l 


88 


F5 89L 




42 92L 


89 


F5 142F 




40 ( )F 


90 


F4 8f 




42 8f 


91 


LO 117L 




32 97L 


92 


4l IF 




hi £ )F 


93 


F5 92L 




42 92L 



By 56, 82 



Test to see if g has converged, 



From 79 
By 57, 83 



From 97 



From 87 



Bare all tne g.*s converged? 
If 80, fcrifcg in print routine. 



From 88, 95 



By 93 



LOCATION 


omm 


N0!EES PAGE 10 


9^ 


F5 ur 








40 IT 




Generate identity aatrix 


95 


LO 117L 








32 92L 




. 


96 


L5 89L 








F^ 117L 




■ 


97 


22 87L 




■ . 


. 


L5 122L 


FrcoR 91 




98 


k2 10^L 








h6 lO^L 






99 


F5 1^2F 








40 2F 






100 


kl 3F 


From Ilk 


- 




F5 2F 




i 


101 


40 2F 
19 ISP 


. 




102 


IA- lCftL 




■ 




k6 lOkh 




, 


103 


L5 30L 
k2 10% 






ICft 


50 ( )F 


By 98, 102 


From 112 




7J ( )f 


%y 9B, 108 




105 


ko F 


1 






50 ( )F 


*jr 103, 109 




106 


75 F 
66 1211* 




- Compute SRG 


107 


S5 F 








>0 (331 )F 


ly 1, 110 




108 


F5 IC&t, 








^2 ICftL 




■ 


109 


F5 10% 
k2 105L 




■ 

^ ■ 


110 


F5 107L 








k2 107L 




■ , 


111 

\ 


F5 3F 







M 10 



LOCJMTXOI 


GREER 


ItfXBB FAfiE 11 




kO 3P 




, ■ . 


112 


LO 2F 








36 1(*L 






113 


IA 2F 
LO U7L 






Ilk 


36 100L 








50 U^L 


; 




115 


26 2QF 




Compute eigeHnrfcliifts and eigenvectors 




26 L 




of ORG 


116 


00 F 








00 SN 


By 51 


'■'•'■ 


117 


80 F 








00 Sk 




118 


OCK-'OO 0032 
0000 0000 J 




119 


80 F 








00 (SK) 


By >9 




120 


oof: 00 ooc& 






6051 7019 J 




121 


OOF 00 0100 
0000 0000 J 




122 


00 1G23S8 






00 s6 






123 


00 t 

00 S9 






124 


92 93^F 




Itekela Qut 10 T*& and stop* -if too 




OF F 




MW^.j^R^erB^are fceing extracted. 




00 15JC 








If 1 8F 








26 2^8F 


. 






26 15H 








00 10K 









09 10F 








L5 128F 


. ■ 





M 10 



Locmosr 


OEDEK 


itoees pact: i? 


i 


00 IP 
kO IF 








2 


00 2F 
iA IF 








3 


kO 3F 
lA 128F 






Printing interlude f 
•CGKcp'ate parameters necessary for 


k 


kO kF. 
50 191F 






Code PI 


5 


JO 7L 










S5 F 






■ ■ . ■ 


6 


kO 8F 
26 999F. 








7 


00 F 

00 1023F 










26 10H 








00 20K 









kl 7F 

92 983F 






1 


92 lk-3F 
kl 8F 






2 


L5 (S8) 
kO 9F 


By 13, 


From 16 


3 


52 1B3 








50 3L 




Print crut G 


Ij- 


26 959F 
50 9F 






5 


75 9F 
i»0 F 






6 


LO 282F 
66 F 






7 


S5 F 








kO (80 )F 


By 1^ 


OoJipute egmacHiml ities and store* 


8 


92 131F 
L5 57L 






9 


50 8F 







M 10 



Loamoi 


omm 




HOTSS 


PAGE 15 




66 9F 








10 


ki lap 

S5 F 








11 


kO (S8) 
19 18P 


By 12 


Ccwpnte qT 




12 


ju4 IJLLt 

k6 12i 








15 


46 2L 
F5 71* 








14 


k2 7L 
F5.7F 








15 


kO 7F 
LO 2?8F 








16 


56 2L 
92 770F 






. 


IT 


92 143F 

L5 (80 )F 


From 22 
By 20 






18 


52 sk 










50 18L 




Erint out cammumlities 


■ 


19 


26 959F 
92 151F 








20 


F5 17L 
k2 17L 








21 


F5 Qf 

40 8F 




. 




22 


LO 2T8F 
52 17L 








25 


92 143F 










F? (551 )F 


By 55, 


Fr$m J> 




24 


56 25L 
26 55L 








25 


L5 25L 
42 26L 


From 24 






26 


42 52L 
L5 ( )F 


By 25 







LOCATION 


OKDER 


j 

SOTSS PAGE 14 


27 


40 4F 








36 28L 




: .'■'■' 


28 


52 2S3 








50 28L 




Brint out largest latent roots 


29 


26 959F 








4l F 






30 


92 131F 

L5 4f 




.. 


31 


LO 282F 








50 31L 


, 




32 


26 152F 




■ •. ■ . 




40 { )F 


By 26 


Compute /X.^1 and store 


33 


F5 23L 
42 23L 


From 24 


: . ■ ■ 


3^ 


F5 lOF 
40 10F 






35 


LO 278F 
32 23L 






36 


4l 8F 
92 l43F 




' 


37 


50 S8 




Constant 




■ L5 37L 


i&m 56 


X. ^ ■ 


36 


46 43L 








F? (331 )F 


By 5^ 




29 


36 40L 
22 53L 


i 


: ■ : . 


40 


L5 38L 
**2 42X* 


From 39 




*a 


Ik 229F 
42 44L 






42 


4i 9F 


■ 






50 { )F . 


From 40, 52 


43 


TJ" ( )F 


By 38, 48 




40 F . 






44 


50 F 
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LOCATION 



ammmaem 



k*j 
kQ 
k-9 
50 
51 
52 
53 
5^ 

I 55 

56 
57 



MMMMMMMMMM 



HOTO 



PAGE 15 



M 10 



75 ( }F 
66 57L 

51 F ; 

52 Sif 
50 ^L 
26 959P 
19 18P 
iA- ^3L 
kQ kjL 
L5 278F 
iA- kkL 
ij-2 1AL 
92 131F 
F5 9F 
ii-0 9F 
LO 278F 
32 42L 
92 135F 
F5 38L 
^2 38L, 
F5 8F ' 
ko 8F 
LO 278F 
32 371* 
OF F 
OOF 00 1000 
0000 0000 J 
26 20H 



From ^1,50 



Cts&renjfce factor loadings! 



4.JL oT x 



Print out factqcr loadings 



From 39 



M>MWM> 




'■ APPROVED BY w 



